Analytic treatment of leading-order parton evolution equations: theory and tests 
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We recently derived an explicit expression for the gluon distribution function G(x, Q 2 ) = xg(x, Q 2 ) 
in terms of the proton structure function F£ f '(x,Q 2 ) in leading-order (LO) QCD by solving the 
LO DGLAP equation for the Q 2 evolution of {x,Q 2 ) analytically, using a differential-equation 
method. We showed that accurate experimental knowledge of F2 P (x,Q 2 ) in a region of Bjorken x 
and virtuality Q 2 is all that is needed to determine the gluon distribution in that region. We re- 
derive and extend the results here using a Laplace-transform technique, and show that the singlet 
quark structure function Fs(x, Q 2 ) can be determined directly in terms of G from the DGLAP gluon 
evolution equation. To illustrate the method and check the consistency of existing LO quark and 
gluon distributions, we used the published values of the LO quark distributions from the CTEQ5L 
and MRST2001LO analyses to form F^ p (x, Q 2 ), and then solved analytically for G(x,Q 2 ). We 
find that the analytic and fitted gluon distributions from MRST2001LO agree well with each other 
for all x and Q 2 , while those from CTEQ5L differ significantly from each other for large x values, 
x > 0.03—0.05, at all Q 2 . We conclude that the published CTEQ5L distributions are incompatible in 
this region. Using a non-singlet evolution equation, we obtain a sensitive test of quark distributions 
which holds in both LO and NLO perturbative QCD. We find in either case that the CTEQ5 quark 
distributions satisfy the tests numerically for small x, but fail the tests for x > 0.03 — 0.05 — their 
use could potentially lead to significant shifts in predictions of quantities sensitive to large x. We 
encountered no problems with the MRST2001LO distributions or later CTEQ distributions. We 
suggest caution in the use of the CTEQ5 distributions. 

PACS numbers: 13.85.Hd,12.38.Bx,12.38.-t,13.60.Hb 
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I. INTRODUCTION 

In a recent paper [l| , we derived an explicit expression 
for the gluon distribution function G(x, Q 2 ) = xg(x, Q 2 ) 
in the proton in terms of the proton structure function 
F2 P (x,Q 2 ) for 7*p scattering. The result was obtained 
in leading-order (LO) QCD by solving the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation Q 
for the Q 2 evolution of F% p (x, Q 2 ) analytically, assuming 
masslcss quarks and using a differential-equation method. 
The method is model-independent and does not require 
individual quark distributions, so accurate experimental 
knowledge of F^ix, Q 2 ) in a region of Bjorken x and 
virtuality Q 2 is all that is needed to determine the gluon 
distribution in that region. 

In the present paper, we extend the method using a 
Laplace-transform technique. We first re-derive the solu- 
tion for G(x,Q 2 ) in terms of 

F2 p (x, Q 2 ) = E e * x Q 2 ) + Q 2 )] W 

i=l 

using our Laplace transform method. The same proce- 
dure can be used to derive G(x, Q 2 ) from the singlet (S) 
quark distribution function 

n f 

F s (x, Q 2 ) = J2 X M x < Q 2 ) + *( x > ^)] ( 2 ) 

and other combinations of the quark distributions. We 
then obtain an analytic solution for Fs(x,Q 2 ) in terms 



of G(x, Q 2 ) by solving the LO gluon evolution equation 
exactly. These are the principal theoretical results in this 
paper. 

To illustrate the method, we have used the LO 
CTEQ5L @ and MRST2001 Q quark distributions to 
construct the proton j*p and singlet structure functions 
as functions of x and Q 2 , solved analytically for the cor- 
responding gluon distributions, and compared the results 
for G(x, Q 2 ) with those given in the fits. We find that the 
analytic and published distributions agree to high accu- 
racy for MRST2001 LO, but, to our surprise, that they 
are incompatible for CTEQ5L, differing significantly at 
large x for all Q 2 . We find the same situation when we 
use the evolution equations for the singlet quark distri- 
bution or for individual quarks to solve for G(x, Q 2 ). We 
have not encountered these difficulties in somewhat less- 
extensive calculations using the more recent CTEQ6L [j| 
distributions, finding good agreement between the ana- 
lytic and published gluon distributions. 

To investigate the inconsistencies further, we have used 
the DGLAP evolution equation for the non-singlet (NS) 
distributions 



Fns 



"f 

E 

i=l 



x [ui(x, Q 2 ) + Ui(x, Q 2 ) 
-di{x,Q 2 )-di{x,Q 2 )}, (3) 



where the sum is over the up- and down-type quarks and 
antiquarks in each generation with both types active. As 
will be shown, Fns can be used in both leading order 



2 



(LO) and next-to-leading order (NLO) to construct sen- 
sitive tests of the quark distributions. 

The tests based on Fms are satisfied for the MRST2001 
LO [H and CTEQ6L [B( quark distributions which also 
satisfied the gluon tests above to high accuracy. However, 
we find that the LO CTEQ5L and the NLO CTEQ5M 
quark distributions in the MS renormalization scheme fail 
to satisfy the NS evolution equation in the same regions 
at large x where the gluons gave difficulty. This shows 
clearly that these published quark distributions are in- 
consistent at large x. 

The errors in the quark distributions and their incom- 
patibilities with G could well affect predictions for exper- 
imental quantities that are sensitive to large x. We con- 
clude that the CTEQ5L distributions, though very con- 
venient because of the analytic parametrizations given 
3, should not be used to make precise predictions at 
large x. 



Since 

F] P (x, Q 2 ) = ^F s (x, Q 2 ) + If ns (x, Q 2 ), (9) 

we should be able to derive G(x, Q 2 ) consistently in LO 
from either of Eqs. ([5]) or ([6]) provided Eq. ([7]) is satisfied. 

We note that Eqs. ^ and ([7]) continue to hold at next- 
to-leading order with modified quark and gluon splitting 
functions in Eqs. ((H) and ©. These are given in the MS 
scheme in Eqs. (4.7a), (4.7b) and (4.7c) of Floratos et al. 

To simplify the equations, we rewrite Eq. ([5]) as 

TT{x, Q 2 ) = x( G(z, Q 2 )K qg (-) ^ , (10) 

J x \ Z / Z 

where we have introduced TT{x,Q 2 ), which is defined 
as 



II. 



PRELIMINARIES 



We first introduce a quantity T 2 {x, Q 2 ), which depends 
only on the proton structure function F2 P (x,Q 2 ), by 



Similarly, we define TT s{x,Q 2 ) m terms of the singlet 
quark distribution Fs(x,Q 2 ) by 



. _ dF^(x,Q 2 ) 



a s 
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dln(Q 2 ) 

4F™(x, 
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F 2 ^(x,Q 2 )ln. 



F^(z,Q 2 ) F^(x,Q 2 ) 



/>W)(i + f)$}- 



(4) 



We define the corresponding quantities Fs(x,Q 2 ) an d 
^Fns(x,Q 2 ) for the singlet and non-singlet quark distri- 
butions obtained by replacing F^ by F$ or F^s- 

The LO DGLAP evolution equations for the 
proton structure functions F^ix.Q 2 ), Fs(x,Q 2 ), and 
Fns(x, Q 2 ) for massless quarks can be written compactly 
in x space as 



Fs(x,Q 2 ) = 



dz 



2nt — x 
J 4tt 

0, 



G(z,Q')K qg (-J^,(5) 
x\ dz 



[ G(z,Q 2 )K 

J X 



(G) 
(7) 



where the sum over quark charges in Eq. (JSJ) includes both 
quarks and antiquarks, and nf in Eq. is the number 
of active quarks, fjvs vanishes only if both members of 
a quark family arc massless or can otherwise be treated 
as active. The LO g — > q splitting function is given by 



K qg (x) = 1 - 2x + 2x 2 . 



(8) 



TT s {x,Q z ) = [2n f ^j T 2 (x,Q 2 ), (12) 

and analogous quantities with appropriate prefactors for 
other structure functions such as > ^2(3) > an< ^ ^2(3) ■ 
We note also that if we wish to evolve a single massless 
quark distribution in LO, we can simply replace F^ p in 
Eq. (0} by Q t = xq t {x,Q 2 ) and TT in Eq. |TT|) by the 
appropriate 

TT^xM 2 ) = (g) _1 ^(.T,Q 2 ), (13) 

so that 



TTi{x,Q 2 )=x G{z,Q z )K qg {-\—, (14) 

J x \Z ' Z 

with i = u,u, d, d, c, c, s, s. 

Since the gluon content of Eq. fTJJ is independent of 
the type of quark, Eq. (fT4|) has the important conse- 
quence that, in LO, 



TT t (x,Q 2 ) 
TT 3 (x,Q 2 ) 



1, 



(15) 



for all x at a fixed virtuality Q 2 , and any choice of i and j. 
As we will soon see, this turns out to be a very sensitive 
test of compatibility of quark distributions with the LO 
QCD evolution equations. 

An alternate form of this relation, which is readily ex- 
tendable to NLO, follows from the evolution equation for 
the non-singlet distribution F N s(x,Q 2 ), TTns{x,Q 2 ) = 
0. Splitting the sums in Eq. ([3]) into sums over up- 



3 



and down- type quarks, we find that the vanishing of 
J r J r N g(x, Q 2 ) implies that 



dtypc(^,Q , 
utypc(^>Q 2 , 



= 1 



(16) 



for all x at any given virtuality Q 2 . This furnishes us 
with an exact numerical test of NLO perturbative QCD. 
We will later see that, while the MRST2001 and CTEQ6 
parton distributions satisfy our constraints for all Q 2 and 
x, the LO CTEQ5L and NLO MS CTEQ5M quark dis- 
tributions satisfy these tests for small x, but fail badly 
at large x. 



III. SOLUTION OF THE LO EVOLUTION 
EQUATION FOR THE PROTON STRUCTURE 
FUNCTION F^ p (x,Q 2 ) 

We now re-derive our analytic solution [l[ of Eq. (JTUJ) 
for the LO gluon distribution, G(x,Q 2 ) using a new 
Laplace-transform method. Introducing the coordinate 
transformation 



v = ln(l/x), 

we define functions G, K qg , and T in v-space by 

G(v,Q 2 ) = G{e-\Q 2 ) 

K qg (v) ee K qg (e~n 
F{v,Q 2 ) = TT{er\Q 2 ). 
Explicitly, from Eq. ([5]), we see that 

K qg {v) = 1 - 2e-' v + 2e~ 2v . 
Further, we can write 



(17) 



(18) 
(19) 



T(v,Q 2 



G(w, Q 2 )e-^ v - w ^ K qg {v ~ w) dw 



G{w,Q 2 )H(v -w)dw, 



where 



H(v) 



e- v K qg {v) 



2e 



2e 



-3v 



(20) 



(21) 



We introduce the notation that the Laplace transform 
of a function H(v) is given by h(s), where 

poo 

h(s) = C[H(v):s}= H(v)e- sv dv, (22) 
Jo 

with the condition H(v) =0 for v < 0. The convolu- 
tion theorem for Laplace transforms relates the transform 
of a convolution of functions G and H to the product of 
their transforms g(s) and h(s), so that 



C 



G{w)H(y — w) dw; 



g(s) x h(s), (23) 



or conversely, the inverse transform of a product to the 
convolution of the original functions, giving 



C-^s) x h(s);v] 



G(w)H(v — w) dw 



H(w)G{v-w)dw. (24) 



Equation (|2"3"|) allows us to write the Laplace transform 
of Eq. (HOI) as 



where 



f( Sl Q 2 ) = g( S ,Q 2 )xh(s), 



f(s,Q 2 )=£[f(v,Q 2 );. 



(25) 



(26) 



and g(s, Q 2 ) and h(s) are the transforms of the functions 
G and H on the right-hand-side of Eq. (|2"0|) . Solving 
Eq. (g5]) for g, we find that 



7(s,Q 2 ) = /( S ,Q 2 )/Ms)- 



(27) 



We will generally not be able to calculate the inverse 
transform of g(s, Q 2 ) explicitly, if only because f(s, Q 2 ) is 
determined by a numerical integral of the experimentally- 
determined function T{v,Q 2 ). However, if we regard 
g(s,Q 2 ) as the product of the two functions /(s, Q 2 )and 
h~ 1 (s) and take the inverse Laplace transform using the 
convolution theorem in the form in Eq. (|24[) and the 
known inverse C~ 1 [f(s,Q 2 );v] — J-(v,Q 2 ), we find that 

G(v,Q 2 ) = C^lfis^^xh-^s^v} 

T(w,Q 2 )J{v~ w)dw (28) 



where J(v) is a new auxiliary function, defined by 

J{v)=C- 1 [br' L {8);v]. (29) 

The calculation of h(s) and the inverse Laplace trans- 
form of of /i _1 (s) are straightforward, and we find that 



J(v) = 36(v) + 5'(v) 
' 6 



-3u/2 



Sill 



+2 cos 



V7 



V7 



(30) 



We therefore obtain an explicit solution for the gluon 
distribution G(v,Q 2 ) in terms of the integral 



G(v,Q 2 ) 



F(w, Q 2 )J(v — w) dw 



= 3T(v,Q 2 



dT{v,Q 2 ) 



dv 



-3(v—w)/2 



4 




dw. (31) 



Transforming Eq. (|3~lj) back into x space, we find that 
the gluon probability distribution G(x, Q 2 ) is given in 
terms of the function TT, assumed to be known, by 



G{x,Q 2 ) = ZTT{x,Q 2 ) 




(32) 



It can be shown without difficulty that this expression 
for G is equivalent to that which we derived in [l| , where 
we first converted the evolution equation for F^ p (x,(5 2 ) 
into a differential equation in v, and then solved that 
equation explicitly. 

Even if we cannot integrate the last term of Eq. |32|) 
analytically, the usual case, it can be integrated numer- 
ically, so that we always have an explicit solution which 
can be evaluated to the numerical accuracy to which 
J-T(x,Q 2 ) is known. We emphasize that this solution 
for G is derived from F^ v on the assumption that the 
quarks are either massless, or effectively so, a situation 
that holds in general for Q 2 ^> AM 2 , where mass effects 
due to the i th (massive) quark arc negligible, or for a 
treatment of mass effects such as that used in CTEQ5 
[H as discussed in Sec. IVl 1 

Leading-order solutions for G(x, Q 2 ) of exactly the 
same form follow from the evolution equations for the 
singlet structure function Fs and the individual quark 
distribution functions, with TTix, Q 2 ) replaced in Eq. 
([52")) by TT s {x, Q 2 ) or TT t {x, Q 2 ), as noted earlier. The 
latter can be used to obtain expressions for G(x, Q 2 ) in 
terms of the structure functions for weak scattering pro- 
cesses, allowing the incorporation of other data sets. 



IV. SOLUTION OF THE LO GLUON 
EVOLUTION EQUATION 

To obtain an expression which relates the singlet struc- 
ture function Fs(x, Q 2 ) indirectly to experiment, we use 
the DGLAP equation for the Q 2 evolution of the gluon 
distribution G{x,Q 2 ). This relates Fs to G which is 



determined by F^ v as shown in the preceding section. 
Again, we simplify the notation by writing this equation 
in terms of the (rescaled) q — > g splitting function K gq (x) 
and a quantity QQ(x, Q 2 ) as 

gg(x,Q 2 ) = x F s (x,Q 2 )K gq (x/z)—. (33) 

J X z 



Here 
and 

gg&Q 2 ) 



K gq (x) 



2 

x 



2 + x, 



(34) 



dG{x,Q 2 ) 
d\nQ 2 



G(x, Q ) [ 3 (33 - 2n/) + 12 ln(l - x) 



Ylx 



x 



G{z,Q 2 ) ~-2 + -- - -7 



z \ z 
z \ dz 



2 \ dz 



+ j\ G (z,Q 2 )-G(x,Q 2 ))(-^-) ^jj, (35) 

for nj effectively massless active quark flavors. As before, 
going to v space, we define the quantities 

F s (v,Q 2 ) ee F s (e- V ,Q 2 ), 

g{v,Q 2 ) = gg{e~\Q 2 ), 

H gq (v,Q 2 ) ee e- v K gg (e- v ) 



e~ v (2e v ~2 + e 

-2v 



= 2-2e~ v 



(36) 



The evolution equation then involves a convolution of Fs 
and H gq . 

Proceeding as before, we take the Laplace transform 
of both sides, solve for the transform fs of Fs, invert the 
transform using the convolution theorem again, and find 
that 



F S (v,Q 2 ) = 



dg(v,Q 2 ) 



dv 



6 

-2 cos 



sin 



V7 



g(w,Q 2 )e 



(y — w) 



-3(v-w)/2 



V7, 

— {v-w) 



dw. 



(37) 



Transforming back into x space, Fs(x, Q 2 ) is given by 

ogg{x,Q 2 ) 



Fs{x,Q' 



dx 

gg{ z ,Q 2 )Q x 



1 We will discuss mass effects in detail elsewhere. 



6 . 
—= sin 
V7 


V7 


2 


—2 cos 


~V7 


2 



In- 



dz 



(38) 
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QQ{x,Q 2 ) is a known function of the gluon distribu- 
tion, G(x, Q 2 ), which, in turn, is a known function of the 
proton structure function F^ix, Q 2 ) that is measured in 
deep inelastic scattering. Fs(x,Q 2 ) is therefore related 
back to K7 p (x, Q 2 ) through G(x, Q 2 ) and the analysis in 
the preceding section, and could be used, in turn, to de- 
termine G(x, Q 2 ) through Eq. |32|) with TTs replacing 
TT. 

It is important to recognize that our method is essen- 
tially model- independent. When different parametriza- 
tions in x and Q 2 are used to fit experimental F^ p data, 
the LO gluon and F$ distributions derived from the dif- 
ferent fits must agree with each other to the accuracy 
of the fits in the regions in which the data exist. The 
results can only depend on the functional forms used to 
fit F2 P (x,Q 2 ) when they are extrapolated to x and Q 2 
outside that region. 



V. NUMERICAL ILLUSTRATION OF THE 
METHOD, AND TESTS OF LO GLUON AND 
QUARK DISTRIBUTIONS 

To test of our methods, we have applied them to pub- 
lished parton distributions, using the quark distributions 
to construct F^ v and Fg and then solving Eq. ([32)) to 
find G. The methods work well, and reproduce the pub- 
lished gluons distributions except in the case of CTEQ5 
Q. In that case, we found to our surprise that there are 
problems at high x with the LO CTEQ5L distributions. 
While these are superseded by more recent CTEQ distri- 
butions 0, 0] , they seem still to be used in some calcula- 
tions, most likely because they have been parametrized 
analytically in a form convenient for calculation. Similar 
tests applied to the MRST2001 Q and the CTEQ6L % 
LO distributions uncovered no problems. 

The results illustrate the sensitivity of the analytic 
methods in testing parton distributions, and may be of 
wider use for that purpose, as well as for deriving G di- 
rectly from data as originally proposed in [l[ . We present 
them here as a demonstration, and as a caution. 

We begin by illustrating the use of the analytic ex- 
pression in Eq. (|3"12|) to derive G(x,Q 2 ) from F^ p (x,Q 2 ) 
in the case of CTEQ5L Q. We take the published LO 
CTEQ5L quark distributions as our basic input, and 
use these distributions to calculate the proton struc- 
ture function F^ p (x, Q 2 ) = £V e 2 x[q t (x, Q 2 ) + q t (x, Q 2 )] 
needed in TT{x, Q 2 ) in Eq. QIJ). Next, we solve Eq. (32|) 
for the G generated by this F^, labelling the solution 
GAnaiyticf^j Q 2 ), and compare the results with the pub- 
lished gluon distributions. 

We follow the procedures of the CTEQ5 group in these 
calculations. The parton splitting functions used are 
those for massless quarks. Quark mass effects are in- 
cluded only in an approximate, x-independent way, with 
xqi and xq~i taken as zero for a massive quark i when 
Q 2 < AI 2 , and the quark treated as fully active when 
Q 2 > Mf, with qi and q~i included in the calculations 



of the structure functions F^ p and Fs- The charge fac- 
tors X)i e ? m Eqs. ([5]) and (TTTj) , and the numbers of 
active quarks nj in Eqs. ©, (fT2")h and later in Eq. (|50")) . 
then change discontinuously at each threshold, but re- 
main constant between thresholds. Despite these discon- 
tinuities, the quark and gluon distributions are continu- 
ous. 2 

Because the mass effects do not depend on x in this ap- 
proach, the evolution equations retain the massless form 
between thresholds, and we can simply use the results de- 
rived above for massless quarks in our analysis, treating 
the distinct inter-threshold regions in Q 2 separately. 3 

In these calculations, we use the LO form of a s (Q 2 ) 
used in CTEQ5L [|, 



a s (Q 2 



-in 



p \n(Q 2 /A 2 )' 
2 

11- 



(39) 
(40) 



with n f = 5 and A 5 = 146 MeV for Q > 4.5 GcV, n/ = 4 
and A 4 = 192 MeV for 1.3 GeV < Q < 4.5 GeV, and 
n f = 3 and A 3 = 221 MeV for Q < 1.3 GcV. These pa- 
rameters give a continuous a s (Q 2 ) with a s (M%) = 0.127 
1- 

The same methods can be used in LO to determine 
the gluon distributions Gi generated by individual quark 
distributions, or by other combinations of the quarks, in 
particular, by combinations of the massless quarks u, d, s 
and by the non-singlet distributions, as will be seen be- 
low. 

Following these procedures, the gluon distribution 
G Analytic^ Q 2 ) we derive analytically from the calcu- 
lated structure functions F^ p or F$ , or from other quark 
combinations, should agree to the accuracy of the calcu- 
lations with GcteQ5l(x,Q 2 ) at all x and Q 2 . 

The results of the calculations based on F^ are shown 
in Fig. [T] The left-hand column compares GAnaiytic (blue 
dashed curves) with Gcteqsl (red solid curves) at Q = 
5, 20, and 100 GeV 2 , (a), (b), and (c), respectively. All 



2 The quark and gluon distributions arc continuous solutions of 
the evolution equations, which are first-order differential equa- 
tions in InQ 2 and integral equations in x. The discontinuities on 
the right-had sides of Eqs. J5]l and JSJl as nf changes at a thresh- 
old are reflected in the final results by discontinuous changes in 
the derivatives of the structure functions dF/dlnQ 2 on the left- 
hand sides of these equations. In particular, the heavy quark 
distributions xqi(x,Q 2 ) and xcji(x,Q 2 ) are identically zero for 
Q 2 < M 2 , then initially rise linearly from zero with increasing 
InQ 2 for Q 2 > M 2 . The contributions of d(xqi)i/d In Q 2 and 
d(xqi)i/dlnQ 2 to dF/dlnQ 2 , identically zero for Q 2 < M 2 , 
therefore jump to non-zero values for Q 2 > M 2 . 

3 Mass effects are treated the same way by CTEQ6 0. The 
more accurate treatments of heavy-quark masses in MRST2001 
[3 and the recent CTEQ analyses beginning with CTEQ6.5 
involve x-dependent effects. These will not affect the tests of 
the MRST2001 distributions given later using only the massless 
quarks. 
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FIG. 1: Left -hand column: comparison of the GAnaiytic 
(blue dashed curves) obtained by solving Eq. (|32[) using the 
CTEQ5L (x,Q 2 ) as input, with Gcteqsl (red curves). 
The u-quark distributions U(x,Q 2 ) = xu(x,Q 2 ) (purple dot- 
dashed curves) are shown to give an independent scale at the 
same x and Q 2 . (a) Q 2 = 5 GeV 2 , n f = 4. (b) Q 2 = 20 
GeV 2 , rif = 4. (c) Q 2 = 100 GeV 2 , n f = 5. Right-hand 
column: the ratio re = GAnaiytic /Gcteqsl at the same Q 2 
as (a), (b), (c). The ratios should equal 1 for all x and Q 2 . 

of the active, physically-relevant quarks are included in 
the input at each value of Q 2 , with rif — 4 for (a) and (b) 
with u, d, s, c active, and n/ = 5 in (c) with b now also 
active. For a scale comparison, we include the CTEQ5L 
up-quark distribution, U(x,Q 2 ) = xu{x,Q 2 ). 

The analytic and fitted G s do not agree well at large x, 
with significant differences between them on the scale of 
G and of the u-quark distribution for 0.05 < x < 0.6. We 
do not believe that the calculations are reliable beyond 
x ?s 0.6: G is very small , but the input u and d distribu- 
tions and the individual integrals in TT are not, and the 
calculation in Eq. (|32|) involves large cancellations and 
becomes sensitive to small errors in the inputs. 

The same discrepancies between the two Gs are il- 
lustrated differently in the right-hand column of Fig. Q] 
where we plot the ratio 

TG = G A nalytic(a;, Q 2 )/G C TEQ5L (x, Q 2 ). (41) 

In all cases, we find that tq is essentially constant and 
equal to one for all x < 0.01 for small virtuality, and for 
all x < 0.05 for large virtuality, showing good agreement 
between analytic and fitted gluon distributions at small 
x. There are again significant deviations from one at 
larger x, contrary to theoretical requirements. 

We emphasize that the reason we find negative values 



of GAnaiytic or re for some x ^ 0.5 is not that -Fr7 p (x, Q 2 ) 
goes negative, but rather that TF{x,Q 2 ) goes negative. 
As can be seen from Eq. (j4|, this is a sensitive function 
of the difference between dF2 P (x,Q 2 )/d\n(Q 2 ) and the 
LO convolution integral of F2 P (x,Q 2 ). In spite of the 
fact that the CTEQ5L quark distributions that go into 
F2 P (x,Q 2 ) are all positive, the resulting combination 
J-J-(x, Q 2 ) becomes negative for large x. The appearance 
of negative values for the gluon distribution function sug- 
gests strongly that the differences between GAnaiytic and 
Gcteqsl in Fig. Q] result from problems with the input 
quark distributions, hence with GAnaiytic, and are less 
likely to arise directly from problems with Gcteqsl- 

We find very similar deviations if we use the singlet 
structure function Fs(x, Q 2 ) in Eq. ([2]) as input, and then 
solve for G(x,Q 2 ) using Eq. ([52")) with TTs replacing 
TT. 

The discrepancies between the analytic and fitted Gs 
are similar in magnitude to the differences between the 
CTEQ5 and more recent CTEQ and MRST gluon dis- 
tributions in this region, and to the changes in the dis- 
tributions that resulted from the addition of inclusive jet 
data to those analyses. Changes of this size are clearly 
significant. 

Since the major contributions to F2 P (x,Q 2 ) at large 
x are from the valence quarks u and d, we have done a 
separate set of calculations using the singlet combination 
Fs,ud = x(u + u + d + d) and Eq. (p{!?|) with TTs.ud re- 
placing TT . No mass effects need to be included in this 
case. The results of the calculations, called G u d(x,Q 2 ), 
are shown in Fig. H for both CTEQ5L and MRST2001 
LO. 

Discrepancies between Gcteqsl and the G analyti- 
cally reconstructed from Fs, u d, here labeled G u d(x, Q 2 ), 
are again evident in the figure. The discrepancies have 
the same pattern and are roughly the same sizes as those 
seen in Fig. [TJ suggesting that the problems originate 
with the u and d distributions that drive the initial 
DGLAP evolution at large x. The results for MRST2001 
LO show agreement between the calculated G u d and 
Gmrst20oi at the level of accuracy of the calculation, 
which required the parametrization of numerical data 
from @. 

To exclude the possibility that the problems with 
CTEQ5 can be eliminated by going to NLO perturba- 
tive QCD, we also looked at the evolution of the non- 
singlet structure function Fms{x,Q 2 ), and at the gluon 
distribution from the evolution of the singlet structure 
function F s (x,Q 2 ), using the CTEQ5M NLO distribu- 
tions as input. We found a similar x-dependence for the 
gluon ratio obtained from the NLO singlet distribution, 
hence, continuing problems. The NS combination, for 
both LO and NLO, is discussed in the next section, and 
shows clearly that there arc problems with the CTEQ5L 
quark distributions. 

We note finally that, to reduce the possibility that 
the discrepancies that we have encountered are a arti- 
fact of the published parametrization of the CTEQ5L 
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FIG. 2: Plots of the n f = 2 G ud (a:,Q 2 ) (blue dashed 
curves), Gfitted(a;, Q 2 ) (red curves), and the n-quark distri- 
bution U(x,Q 2 ) = xu(x,Q 2 ) (purple dot dashed curves) vs. 
x, for x > 0.1. Left-hand column: Gfittcd = Gcteqsl distri- 
butions for (a) Q 2 = 5 GeV 2 , (b) Q 2 = 20 GeV 2 , (c) Q 2 = 100 
GeV 2 . Right-hand column: Gflttcd = Gmrst20oi LO distribu- 
tions for (d) Q 2 = 5 GeV 2 , (e) Q 2 = 20 GeV 2 , 
GeV 2 . 



(f) Q 



100 



parton distributions Q , we have also compared them to 
the corresponding numerical CTEQ5L distributions from 
the Durham web site Q and concluded that they are nu- 
merically compatible. 

Because the GAnaiytic^, Q 2 ) derived from either 
F^ p (x,Q 2 ) or F s (x,Q 2 ) fails to match C C teq 5 l(x, Q 2 ) 
well at large x, we cannot expect the inverse problem, 
the analytic derivation of Fs from G to work well for the 
CTEQ5L distributions, and do not include such calcula- 
tions here. 



TESTS OF LO AND NLO NON-SINGLET 
DISTRIBUTIONS 



VI. 



For massless quarks in LO, the relation in Eq. fT5 
holds for any quark or anti-quark pair i and j, 



TT % (x,Q 2 ) 
TTj{x,Q 2 ) 



(42) 



As discussed earlier, this can be generalized read- 
ily to NLO using the non-singlet evolution equation 
Fns(x, Q 2 ) = 0, Eq. 0with the NLO NS splitting func- 
tions of Floratos et al. [6j . If we separate the sum in Eq. 



Q into sums over u- and d-type quarks, defined as 

u typo = x[u(x, Q 2 )+u(x, Q 2 ) +c(x, Q 2 ) +c(x, Q 2 )] , (43) 
d type = x[d(x, Q 2 )+d(x, Q 2 ) +s(x, Q 2 ) +s(x, Q 2 )] , (44) 

this relation becomes 



TT N s{x 7 Q 2 ) = FF ut y pc (x, Q 2 ) - TT dtypc {x, Q 2 ) = 0, , 

(45) 

where the two terms are to be calculated separately using 
the same NLO splitting functions as for the full Fns- 
Dividing by the it-type term, we obtain the ratio 



J~ $~ dtype(^, 



utypo 



(x,Q 2 ) 



= 1: 



(46) 



a relation true for both LO and NLO when all the rele- 
vant quarks are active. 

To use this relation to check the consistency of the 
quark distributions, we insert the appropriate LO or NLO 
quark and anti-quark distributions in Eqs. (|4"3"|) and(T4"4"]) 
and use them to generate rNs(a;, Q 2 ) numerically using 
the LO or NLO splitting functions. 

In the case of the NLO CTEQ5M MS distributions, we 
use 



a s {Q 2 ) 



4tt 



p \n(Q 2 /A 2 ) 

2ft ln[ln(Q 2 /A 2 )] 



(3 2 ln(QVA 2 ) 



ii - 5 », 

19 

51 ~ 



(47) 
(48) 
(49) 



with a s (M 2 ) — 0.118 for rif — 5, matching to the a's for 
n f = 4 and 3 at Q = 4.5 and 1.3 GeV. 

In the left-hand column of Fig. [3l we plot tns vs. x, for 
Q 2 = 5, 20, and 100 GeV 2 and n f = 4. The solid (red) 
curves were calculated using the LO CTEQ5L quark dis- 
tributions as input, and the dashed (blue) curves, the 
NLO CTEQ5M quark distributions in the MS renormal- 
ization scheme. We see that the patterns are very similar 
for all Q 2 , with ratios tns ~ 1 in absolute normalization 
for all x < 0.01, and with very large deviations occurring 
at large x. This result, independent of the value of Q 2 , 
shows explicitly that there are problems with the quark 
distributions at large x: the NS evolution equation is not 
satisfied in either leading or next-to-leading order. 

To look at the size of the deviations in a more quanti- 
tative way, we introduce the notion of a LO "non-singlet 
gluon" distribution Gns(x,Q 2 ), defined by the equation 

F NS (x,Q 2 ) = 2n f ^-x J G NS (z,Q 2 )K qg {-) — . (50) 

Gns = if !Fns(x,Q 2 ) is a solution of the non-singlet 
DGLAP equation for an even number of effectively mass- 
less quarks. Further, Gns is determined mainly by the u 
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vanish indicates an inconsistency of the u and d distri- 
butions as solutions of the DGLAP evolution equations, 
and absolute magnitudes of the ratio near 1 would indi- 
cate that the integrated discrepancies are large. 

As seen from the right-hand column in Fig. [3J the ratio 
r NS,G is zero for small x, but has large absolute values 
for x > 0.1 for Q 2 = 5, 20, and 100 GcV 2 , Figs. [31(d), 
(c), and (f) respectively. In particular, around x ~ 0.3, 
the ratio is about unity, indicating a serious discrepancy 
near the maxima of the valence-quark distributions. 

The pattern in x of the discrepancies between the gluon 
distributions found analytically from either or F$ 
and the distributions published by the CTEQ5 group 
closely follows the pattern seen in the departure oITTms 
from zero. For example, TTms and the differences be- 
tween the analytic and fitted gluon distributions all be- 
come negative or positive in the same regions. This would 
seem to indicate that the quark distributions at large x 
(and, in particular, the dominant u and d valence distri- 
butions) are the origin of the difficulty, the combination 
of their shapes and Q 2 dependence not being compatible 
with the DGLAP evolution equations, either in LO or 
in NLO in the MS renormalization scheme. 



FIG. 3: Left-hand column: the LO and NLO non-singlet 
d-type to it-type ratio of Eq. (|46|l for nj = 4 and Q 2 = 5, 
20, and 100 GeV 2 , (a), (b), and (c), respectively. The ratio 
should be 1 at all x and Q 2 . The solid (red) curves are LO, cal- 
culated from Eq. (@3j), Eq. (g4j), and Eq. Jl6} using CTEQ5L 
data for quark distributions . The dashed (blue) curves are 
NLO calculated for the MS CTEQ5M quark distributions 0|. 
Right-hand column: rjvs.G = C7jv S (a;)/GcTEQ5L(a;) vs. x, for 
Q 2 = 5, 20, and 100 GeV 2 , (d), (e), and (f), respectively. 
The non-singlet "gluon" distribution Gns {x) calculated from 
u and d quarks is defined in the text in Eq. (|50|l and is ex- 
pected to be zero. 



and d distributions at large x, so measures the accuracy 
of those distributions. 

The normalization on the right-hand side of Eq. ([501) 
has been chosen so that this equation is the formal coun- 
terpart of the singlet equation Eq. (JJ). 

The ratio 



VII. 



CONCLUSIONS 



TNS.G 



G NS (x,Q 2 ) 

G(x,Q 2 ) ' 



(51) 



should vanish for solutions of the DGLAP equations, and 
any deviations are scaled to the size of the actual (singlet) 
gluon distribution Gcteqsl- 

We consider specifically the non-singlet combination of 
massless u and d quarks, 

F N s,ud( x iQ 2 ) = x[u{x,Q 2 ) +u(x 1 Q 2 )] 

-d{x,Q 2 )-d{x,Q 2 )], (52) 

and calculate r/vs.G for that combination, shown in the 
right-hand column of Fig. [3J Any failure of the ratio to 



We have first developed a powerful method for 
the analytic solution of the LO DGLAP evolu- 
tion equations based on Laplace transforms. This 
method allows us to determine G(x,Q 2 ) directly 
from F2 P (x,Q 2 ) or other measurable structure 



functions such as F, 



jZ 



2(3) 



F z 



and Fff a \ , provid- 



ing close and simple connections to experiment for 
massless or effectively massless quarks. 

We can also determine G(x, Q 2 ) analytically from 
the singlet quark distribution F${x,Q 2 ) if the 
quark distributions are known, or Fs{x,Q 2 ) from 
G(x, Q 2 ) when the latter is known. The set of rela- 
tions provide consistency checks on LO quark and 
gluon distributions obtained in other ways. These 
are the principal theoretical results of the paper. 

As an illustration of our methods, we have used the 
analytic solutions to the evolution equations to ob- 
tain tests of the consistency of published quark and 
gluon distributions. In particular, we compare the 
gluon distributions determined analytically from 
structure functions F2 P (x,Q 2 ) and Fs(x,Q 2 ) cal- 
culated from the quark distributions determined in 
various analyses, to the gluon distributions given 
by the same analyses. We have found no prob- 
lems for analytic and fitted gluon distributions for 
the MRST2001 LO or CTEQ6L solutions, while the 
analytic and fitted gluon distributions for CTEQ5, 
though consistent at small x, show small, but sig- 
nificant, deviations in the large x region, indicat- 
ing that the quark and gluons distributions are not 
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completely consistent (lOj. 

3. A further analysis using the non-singlet structure 
function F NS (x,Q 2 ) and both the LO CTEQ5L 
and the NLO CTEQ5M quark distributions showed 
definitively that there are problems with those dis- 
tributions at high x. The quark distributions are 
not consistent with either LO (CTEQ5L) or NLO 
perturbative QCD in the MS scheme (CTEQ5M) 
in the sense that they do not satisfy the appropri- 
ate non-singlet relations for x >, 0.05 even though 
they satisfy them very well at small x. The dis- 
crepancies, of unknown origin, are large for x ~ 0.3, 
and could have serious consequences for predictions 
of processes sensitive to that x region. We con- 
clude that the CTEQ5 distributions 0] should not 
be used to make such predictions. Again, we found 
no problems with the MRST2001 distributions 
or the more recent CTEQ6 distributions [1, @| ■ 



4. We suggest that groups that calculate parton dis- 
tributions numerically from the coupled DGLAP 
equations construct the appropriate analytic solu- 
tions to test the consistency of their results. 

We note, finally, that we are working on an extension of 
our analytic LO solution for G(x,Q 2 )to include massive 
c and b quarks, using the methods of [7|, and on the NLO 
gluon solution. 
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